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A point particle approximation to the classical dynamics of well separated vortices of the 
abelian Higgs model is developed. A static vortex is asymptotically identical to a solution of the 
^vq linearized field theory (a Klein-Gordon/Proca theory) in the presence of a singular point source 

at the vortex centre. It is shown that this source is a composite scalar monopole and magnetic 
£N| , dipole, and the respective charges are determined numerically for various values of the coupling 

^ ■ constant. The interaction potential of two well separated vortices is computed by calculating 

' the interaction Lagrangian of two such point sources in the linear theory. The potential is used 

. to model type II vortex scattering. 

m 

1 Introduction 

The abelian Higgs model |jj is a relativistic field theory consisting of a complex scalar field </> 
coupled to a U(l) gauge field A^, and given a Higgs symmetry-breaking self interaction which 
allows topologically stable solitons to exist. The present work concerns the (2 + l)-dimensional 
model, in which the solitons are simple lumps of energy called vortices (in (3 + 1) dimensions the 
model admits extended string-like solitons Q whose gravitational effects may be important in early 
Cosmology). Multivortex dynamics falls into one of three regimes, depending on the Higgs mass /i: 
if /i is small, vortices attract one another (the type I regime), while if \x is large they repel (type II), 
■ and at one critical value of /i static vortices exert no net force on one another. 

In this paper, we calculate the interaction potential of two widely separated vortices by means 
of a novel approximation. The idea is that, viewed from afar, a static vortex looks like a solution 
of a linear field theory in the presence of a singular point source at the vortex centre. As will be 
shown, the appropriate point source is a composite scalar monopole and magnetic dipole in a Klein- 
Gordon/Proca theory. If physics is to be model independent, then the forces between well separated 
vortices should approach those between the corresponding point particles in the linear theory as 
the separation grows. Proceeding on this assumption, we calculate the asymptotic static two vortex 
potential. Our answer agrees with that of Bettencourt and Rivers, which was derived using a field 
superposition ansatz ||. 

The monopole charge and dipole moment of the composite point source depend on the Higgs 
mass /i. To fix their values, one must solve the static nonlinear field equations with vortex boundary 
conditions, and such solution is perforce numerical. Using a simple Runge-Kutta scheme we have 
determined these charges for various values of /i 2 . Here our results disagree with the work of 
Bettencourt and Rivers, who, while leaving undetermined the constants analogous to these charges, 
make certain assumptions about them which appear ill justified. 

The asymptotic potential reproduces the aforementioned dynamical trichotomy into type I, 
type II and critical regimes found in the nonlinear model. As an application, the scattering of 
type II vortices is calculated and compared with numerical simulations. 
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2 The abelian Higgs model 



We begin by reviewing some standard results concerning the abelian Higgs model |jj . The Lagrangian 
density is 

c = \d^dh- \f, v f^ - ^(H 2 - i) 2 (i) 

where D^cj) = {da + iA tl )<f> is the gauge covariant derivative, F uv = d^A v — d v A u is the field strength 
tensor, and M 2+1 has signature (+,—,—). Note that the electric charge and vacuum magnitude of 
the Higgs field have been normalized to unity, leaving only one parameter /z, the Higgs mass. With 
these conventions the model is critically coupled if fj, = 1, and in the type I (type II) regime if /i < 1 
(/i > 1). The Eulcr-Lagrange equations derived from C are 

zV^-i(KH 2 -i) = o 

5^ + ^(^0-^0) + |0| 2 ^ = 0, (2) 

a set of coupled, nonlinear, hyperbolic partial differential equations of which no nontrivial solutions 
are known. 

If a configuration is to have finite energy, the fields should satisfy the following boundary condi- 
tions as r = |x| — > oo, 

M-l, D0^O, (3) 

whence </>oo := linv^oo </> takes values on the unit circle in C and is thus a continuous map — > S , 
where represents the circle at spatial infinity. Such configurations fall into disjoint homotopy 
classes labelled by the degree of 4>oo (an integer, n, also called the winding number of <j>). The 
time evolution defined by the field equations conserves energy, and so cannot take a field in one 
homotopy class into a different class, since such an evolution, being continuous, would define a 
homotopy between distinct classes, a contradiction. Hence n is a topologically conserved quantity. 

Another consequence of the boundary conditions is that the total magnetic flux is topologically 
quantized, for 

-J d 2 xF 12 = 2mr. (4) 

follows from (||) and Stokes' theorem. A static solution with n = 1 is a stable lump of energy called 
a vortex. This may be visualized as a single flux tube in R 3+1 with translation symmetry along the 
x axis. It has a total flux of 2ir penetrating the physical plane. 



3 Vortex asymptotics 



The first task in the point particle approximation is to find out what a static vortex looks like far 
from its core || ||. We place the vortex at the origin, and use plane polar coordinates. Substituting 
the ansatz 

= a(r)e lf> 

(A ,A r ,A e ) = (0,0, -a(r)) (5) 
(a is real) the field equations (||) reduce to two coupled nonlinear ordinary differential equations, 
d 2 a 1 da 1 .„ , 1 9 . , . 

S-^ + a-V 2 - 0, (6) 
or r dr 
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the equations for Aq and A r being trivially satisfied. Regularity demands that cr(0) = a(0) = 
while the boundary conditions (|^) become 

lim a(r) = lim a(r) = 1. (7) 

r — >oo r — >oo 

Note that the ansatz has unit winding by construction. No exact solutions of (^|) with these boundary 
conditions are known, but numerical solutions suggest that both a are a are monotonic functions of 
r, and that the ansatz produces an isolated lump like structure. 

We are interested in the asymptotic forms of a and a, and for these explicit expressions do exist. 
Define the functions a and (3 such that 

a(r)=l + a(r), a(r) = 1 + f3(r). (8) 

Then (fjj) implies that a and (3 are small at large r, so we substitute (||) into (|^) and linearize in a 

and f3, 

2 f d 2 a I da \ ^ 
^ \d(fir) 2 fir d(fxr) ) 

dr 2 (r)^~rG?r(r) ( ^ r 2 ) r ^ ^ 

These are the modified Bessel's equations of zeroth order for a in fir and first order for [3/r in r 
respectively. Hence, at large r, 

(3 ~ (10) 

where K n is the n-th modified Bessel's function of the second kind Q. Note that Ki = —K' Q . 

Since we have linearized the field equations, the asymptotic solutions contain unknown scale 
constants q and m which can only be fixed by solving (Q) numerically. Rather than solving the 
boundary value problem c(0) = a(0) = 0, er(oo) = a(oo) = 1, we solve the initial value problem 
<r(0) = a(0) = using ct'(0) and a'(0) as shooting parameters. In fact, this is a slight oversimpli- 
fication: due to the singularities of equations (|^) at the origin, we must shoot from r = ro, some 
small positive number, rather than r — 0. Substituting Taylor expansions for a and a into (^) we 
find that, near the origin, 

a = air+ifli (b 2 + r 3 + 0{r 5 ) 

a = b 2 r 2 -4r 4 + 0(r 5 ). (11) 
o 

We use a\ and hi as shooting parameters, adjusting them until the numerical solution has cr^oo) w 
1 ps a(r oc ), where Too is some large positive number, the effective infinity. Having generated such a 
numerical solution, we compare it at large r to the asymptotic forms (|Io| ) and deduce q and m. 

The results of this procedure using a fourth order Runge-Kutta method with ro = 10~ 8 and 
roc = 10 for various values of fi 2 are presented in table 1. That r^ is so small is unfortunate but 
necessary: at large r the field equations reduce to Bessel's equations, which have two independent 
solutions, one exponentially decaying and the other exponentially growing. We seek to pick out the 
former and completely exclude the latter, an impossible task. Hence, all numerical solutions blow up 
at large r, and even though ai and b 2 were tuned to six decimal places, the Runge-Kutta algorithm 
could not shoot beyond r = 10. 
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M 


q 


TO 


r c 


0.4 


-7.54 


-14.92 


4.23 


0.6 


-8.71 


-12.61 


3.75 


0.8 


-9.70 


-11.31 


3.43 


0.9 


-10.14 


-10.89 


3.22 


1.0 


-10.58 


-10.57 




1.1 


-10.98 


-10.31 


2.98 


1.2 


-11.43 


-10.06 


3.07 


1.3 


-11.80 


-9.85 


2.96 


1.4 


-12.23 


-9.66 


2.95 


1.6 


-13.04 


-9.34 


2.88 


1.8 


-13.97 


-9.09 


2.87 


2.0 


-14.50 


-8.86 


2.72 



Tabic 1: Numerical values of vortex scalar charge q and magnetic dipole moment to. The other data are 
the critical points of the static intervortex potential. 



Nevertheless, the qualitative nature of the /_i 2 dependence of q and m is clear. In particular, two 
points about the numerical charges are noteworthy. First, at critical coupling (/i 2 = 1), q i=s m. In 
fact, one can prove that q = m exactly in this case, because the static /i 2 = 1 vortex satisfies a pair 
of first order field equations. These are deduced by means of an argument due to Bogomol'nyi jfj], 

{D l +iD 2 )<f ) = (12) 
-F l2 + \m 2 -l) = (13) 
(in the A = gauge) and within our ansatz take the form 

r^-(l-oW = (14) 
dr 

-^ + (a 2 -l) = 0. (15) 
r dr 

Substituting (|j) into (|l4| ) yields, on linearizing, 

P=-r%- (16) 
dr 

and so a = qK (r)/2ir (3 = qrKi(r)/2it. Thus m = q. We emphasize that this argument works 
only at critical coupling. 

Second, \q\ and \m\ are monotonic functions of /i 2 , \q\ increasing and \m\ decreasing. Bettencourt 
and Rivers || also find the asymptotic forms (|To|) , but leave their charges analogous to q and to 
undetermined. For purposes of calculation, they make two assumptions about the charges which, 
in the light of table 1, may prove ill-justified. First, they assume that q = to is approximately 
true away from /Lt 2 = 1, whereas in our results, q/m varies between 0.50 and 1.64. Second, they 
impose the condition that the magnetic flux of a vortex should vanish at r — and deduce that 
to — —2~k (it is unclear why a condition at r = should directly constrain the asymptotic behaviour 
as r — > oo). This result may be valid for very large /i 2 , but is certainly flawed close to /i 2 = 1 
since the /i 2 = 1 vortex has maximum magnetic flux at the origin, as is easily seen from the lower 



Bogomol'nyi equation (13). So, they combine assumptions which are individually true only in widely 



disparate physical regimes. 
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4 The point vortex 



The next task is to replicate the vortex asymptotics, found above, in the linear field theory by cou- 
pling the fields in standard fashion to a scalar density p and a vector current j„, as yet undetermined. 
To linearize the abelian Higgs model, we choose gauge so that <f> is real. Defining the field ip = l — (f>, 
the vacuum is then ijj = 0, and the linear Lagrangian density is obtained by expanding (|l|) up to 
quadratic order in ip and A^, 

Arcc = \d^3^ - ~ M V - \f^F^ + \A»A». (17) 
Including the external source Lagrangian density, 

C SOUTce = (ytp-j ll A fi , (18) 
we obtain the following massive, inhomogeneous wave equations for ip and A^, 

(□ + M 2 >0 = p (19) 

All gauge freedom has been exhausted, and there is no global U(l) symmetry of £f r0 o with whose 
Noether current we can identify because tp is real. Hence there is no reason to assume that 
is a conserved current, and we cannot set the extra "fictitious current" term in the Proca equation 
@ to zero. 

To make comparison with the asymptotic vortex fields, these must first be converted to the real 
4> gauge. Since <f> has non-zero winding, there is no gauge transformation regular on all M 2 which will 
accomplish this. However, we only require comparison at large r, so for our purposes it is sufficient 
that the transformation be regular on R 2 \{0}. Since a singular point source will be introduced into 
the linear theory, this is from the outset regular only on R 2 \{0}. So, we unwind the static vortex 
(||) with gauge transformation i— > e~ ,e </>, A^ i— > A^ + 8^8 to obtain 

cf, = a ( r )~ 1 + ±K (nr) (21) 
Ag = -a(r) + l~^rifi(r), (22) 

while A r = Aq = 0. It is convenient to introduce a unit vector k in a fictitious third direction 
perpendicular to the physical plane, so that the R 3 vector product can be defined. In terms of the 
2-vector field A, the unwound asymptotic behaviour is 

Tn vn " 

A~-—K^r)0 = -—kxVKo(r). (23) 

We thus seek sources p and j M such that the solutions of (jl^,|2^) are 

V. = ^-K Q (pr) (24) 
(A°,A) = (0,-^kxVJr o (r)). (25) 
The static Klein-Gordon equation in (2 + 1) dimensions has Green's function Kg, 

(-A + p, 2 )K {pr) = 2tt5(x). (26) 
Substituting (Eij) into ( |l~9| ) and using (||^) one finds that 

p = ? 5(x). (27) 
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Similarly, substitution of ( p5| ) into ( |20| ) yields 



j - V(V • j) = -mk x V5(x). (28) 

Taking the divergence of ( p8[ ) one sees that V • j is a solution of the homogeneous static Klein-Gordon 
equation, so if j is a point source (meaning j = except at x = 0) then V • j = everywhere. Thus 
the unique point source satisfying ( p8| ) is 

j = -mk x Vtf(x). (29) 

Since A = we take j° = 0. The physical interpretation of these expressions for p and j is that 
the point source consists of a scalar monopole of charge q and a magnetic dipole of moment m 
perpendicular to the physical plane. Both q and m are negative (see table 1). We refer to this 
composite point source as the point vortex. 



5 The static intervortex potential 

Having found the scalar charge and magnetic dipole moment carried by a point vortex, it is straight- 
forward to calculate the force between two such vortices held at rest, in the framework of the linear 
theory. The interaction Lagrangian for two arbitrary (possibly time dependent) sources (pi, jm) 
and (/02,i(2)) is 

L int = + L A = J d 2 x Pl 4>2 - J d 2 xf {1) A^ (30) 



where (ipi > ^(i) ) are the fields induced by source (pi , ) according to the wave equations (|19|]2C|) . 
This is found by extracting the cross terms in J d 2 x (£f roc + £ SOU rce) where (p, j) is the superposition 
of the two sources, and (iff, A) is a superposition of the induced fields. The expression (^) looks 
asymmetric under interchange of sources 1 «-> 2, but in fact L lnt is symmetric as may be shown using 
the wave equations ( |l9| , p0| ) and integration by parts. 

Now consider the case of two static point vortices, vortex 1 at y and vortex 2 at z. Then 
pi = qS(x — y), while the scalar field due to p 2 is tp2 = QK (p\x — z\)/2tt. Hence, 

= / d 2 x |^(x - y)K Q (p\ X - z|) = 7^K (p\y - z|). (31) 

The magnetic interaction is similar: jfL = 0, j^) = — mk x V<5(x — y) while A9 2 , = 0, A( 2 ) = 
— mk x VLT (|x - z|), so 

d 2 x— [kx V<5(x-y)] • [kx VJT (|x-z|)] 
Ztt 

m2 a T ^ n in 

= -^A^ (|y-z|) 

2 

= -^o(|y-z|) (32) 



using (26) with y ^ z. The total interaction Lagrangian is a function of |y — z| only, so we interpret 
—Ljnt as the potential energy of the interaction, 

U=^[m 2 K (r)-q 2 K (p J r)] (33) 

where r is the vortex separation, that is r = r(cos$,sini9) := y — z. This is the same potential as 
found in 101, but we arrived at it via a different route. 
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This potential is consistent with the partition into type I, critical and type II regimes. The 
central force due to U is 

- U\r) = -^[m 2 ifi(r) - W 2 ^i(,ir)]. (34) 

If fj, < 1, then K\ — > at large r faster than Ki(fir), so scalar attraction dominates over magnetic 
repulsion and the force is negative, consistent with type I behaviour. If /i > 1, the reverse is true and 
the force is positive at large r, consistent with type II behaviour. Potentials for fi 2 = 0.4 (type I) and 
fi 2 = 2.0 (type II) are plotted in figure 1. At /i = 1, m = q, as explained in section || so U = and 
there is no net force at all. This consistency at large r emerges regardless of the specific values of m 
and q away from ji = 1, and may be attributed to the inverse relationship between a field's mass and 
its range. At moderate r, the /i dependance of q/m becomes important. Given that K\ is a strictly 
decreasing function, it is clear from ( |3~i| ) that there exists a unique critical point of U for each /i 7^ 1 
if and only if m/q > yfyL when fi < 1 and m/q < JJL when /it > 1. Our numerical work suggests that 
m/q easily passes these criteria. The rightmost column of table 1 presents the approximate critical 
vortex separation r c for each value of /i 2 . No such critical points were found by Rebbi and Jacobs 
H , who obtained approximate static intervortex potentials by numerically minimizing the potential 
energy functional subject to the constraint that <j> has two simple zeros separated by a given distance. 
So these equilibria are probably artifacts of our approximation, which we take to break down for 
r < r c . Of course, some kind of breakdown is to be expected: vortices are not point particles, as in 
our picture, and when they approach one another closely enough their overlap produces significant 
effects. 
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Figure 1: The potential function U(r) for fi 2 — 0.4 and /1 2 = 2.0. 



6 Type II vortex scattering 

The interaction potential U provides a very simple model of two-vortex dynamics: the dynamics of 
two point particles each of mass M (the energy of a single vortex at rest, a /z 2 dependent quantity) 
interacting via the potential (|3^). Ignoring the (trivial) centre of mass motion, the Lagrangian of 
such a mechanical system is 

L = ^M{r 2 + r 2 i} 2 ) - U(r), (35) 
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since the reduced mass of the system is M/2, This is a manifestly bad model if /i = 1, because it 
would predict that there is no scattering at all, in conflict with the results of numerical simulations 
H and the geodesic approximation |?]]. Away from critical coupling, one might expect the potential 
U to dominate over velocity dependent corrections, at least at moderately low speeds, so the above 
model, although simple, may give a good quantitative account of long range vortex interactions. We 
choose to study type II vortices because these provide a simple, clear-cut dynamic problem: vortex 
scattering. Type I dynamics is slightly more complicated in that vortices can scatter or form bound 
states depending on the initial conditions. The coupling chosen for the type II numerical simulations 
of U is [i 2 = 2, a choice which we follow for purposes of comparison. In the Lagrangian ([35]), the 
constants q and m are already known for \j? = 2, but the vortex mass M is not. Rather than 
attempt to calculate M from our numerical solution, we use the careful numerical analysis of Rebbi 
and Jacobs Unfortunately, they found M for each of a regular sequence of [i values, rather than 

values, so the /i = \[2 value is not quoted. However, a graph of /i against M is very nearly linear, 
so we use linear interpolation to estimate the fi = 1.41421 . . . mass from the y, = 1.4 and fi = 1.5 
masses given. The result is M — 1.512307T. The potential for this coupling is plotted in hgure 1. 

From the plot of U (r) we see that all trajectories which do not encroach on the interior region 
r < r c are scattering trajectories. By time-translation and rotational symmetries, we can, without 
loss of generality, take the point of closest approach (at which r — tq say) to lie on the d = ray 
and occur at time t = 0. It is then straightforward to show that lim^oo 1? is 



± m r )-U { r)) + J>(^-^ 



(36) 



where J = r 2r d is the conserved angular momentum conjugate to d. The deflection angle 9 is 
7T — 2^00. Solving the scattering problem then amounts to numerically approximating this integral. 
Note that there is an integrable singularity in the integrand at r — r . This presents no problem in 
principle, but it must be treated carefully in any numerical algorithm. Schematically, we handle the 
integral as follows, 





f + 




/ + 


Jr +S 


J A 



« 0* + #nc + 0a, (37) 

where 5 is small (S = 0.1) and A is large (A = 15). The contribution $g is calculated by Taylor 
expansion of the integrand about r = Tq, while i) mm is evaluated using the Newton-Cotes rule. At 
large r the potential falls off exponentially, so for r > A we set {7=0 and calculate $a in the free 
vortex approximation. 

To make comparison with the numerical simulations described in |^] we calculate Q as a function 
of impact parameter b for scattering at impact speeds = 0.1,0.2,0.3 and 0.4. The connexion 
between (6, Woo) and (ro, J), the parameters used in ( |36"| ) is found using energy and angular mo- 
mentum conservation. One might worry that at high Voo and low b the vortices will penetrate the 
r < r c zone and become unrealistically captured. In fact even in a head on collision, the speed 
required for this is greater than 0.4, so the problem is never encountered. The results are shown in 
figure 2. As one would expect, the approximation fares reasonably well for large impact parameters 
and moderate speeds, but less well in scattering processes where the vortices approach one another 
closely. The fit to the numerical simulations of j^] could be improved by adjusting the values of q 
and m, a procedure which we eschew on the grounds that it would corrupt the deductive nature of 
the model. 



8 





01234 01234 
b b 



Figure 2: The scattering of fi 2 = 2 (type II) vortices: deflection angle versus impact parameter b 
at four different impact speeds. The solid curves were produced using the point source approximation, 
the crosses by numerical simulation of the full field equations /p]/. 



7 Conclusion 

In this paper we have presented a point source formalism for long range vortex dynamics. We used 
this framework to rederive the static intervortex potential from a new perspective and solved the 
scattering problem for fi 2 = 2 vortices, finding reasonable agreement with numerical simulations, 
despite the simplicity of the mechanical model. It would be straightforward to apply the method 
to other situations of interest: to derive the asymptotic forces between a vortex- antivortex pair, or 
higher winding conglomerations (in the type I regime), or larger collections of vortices for example. 

A less straightforward extension of the present work concerns the long range interactions of 
critically coupled vortices. These exert no forces on one another when at rest, but do affect one 
another when in relative motion. So the scattering of critically coupled vortices is highly nontrivial, 
as evidenced by numerical simulations || and a partly numerical implementation of the geodesic 
approximation §. In the case of Yang-Mills-Higgs theory, another field theory the scattering of 
whose critically coupled solitons (BPS monopoles) has been extensively studied using the geodesic 
approximation, Manton has devised a method for finding long range velocity dependent forces within 
the point particle approximation [fl0[ . The idea is to calculate the interaction of one point source with 
the retarded potential generated by another moving along some trajectory. For BPS monopoles, this 
led to a formula for the asymptotic metric on the two monopole moduli space, which turned out to 
be in precise agreement with the asymptotic form of the Atiyah-Hitchin metric . The method has 



been adapted to several models |12, [13|, but always where the linearized theory is massless. In the 
vortex case, the linearized theory is massive, so the nontrivial problem is to find a suitable substitute 
for ordinary retarded potentials (field disturbances no longer travel uniformly at the speed of light, 
so standard retarded potentials are not appropriate). If velocity dependent intervortex forces can 
be derived by this means, one could deduce the asymptotic form of the metric on the two vortex 
moduli space, at present known only numerically. 
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